Phonon-induced dephasing of singlet-triplet superpositions in double quantum dots 

without spin-orbit coupling 
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We show that singlet-triplet superpositions of two-electron spin states in a double quantum dot 
undergo a phonon-induced pure dephasing which relies only on the tunnel coupling between the dots 
and on the Pauli exclusion principle. As such, this dephasing process is independent of spin-orbit 
coupling or hyperfine interactions. The physical mechanism behind the dephasing is elastic phonon 
scattering, which persists to much lower temperatures than real phonon-induced transitions. Quan- 
titative calculations performed for a lateral GaAs/AlGaAs gate-defined double quantum dot yield 
micro-second dephasing times at sub-Kelvin temperatures, which is consistent with experimental 
observations. 

PACS numbers: 73.21. La, 03.65.Yz, 72.10.Di, 03.67.Lx 



i 

G 

O 

o 



> 
m 

00 

p 

cn 

o 

o> 
o 



X 



The idea of encoding quantum information in spin 
states of electrons confined in quantum dots (QDs) is 
considered to be one of the most promising approaches to 
solid-state quantum computing. The original proposal of 
using single spin as quantum bits [l[ was followed by more 
sophisticated concepts in which a logical qubit was to be 
coded in a system of a few spins @, Q . A class of proposed 
implementations uses two-spin states in double quantum 
dots (DQDs), assigning the quantum logical values to sin- 
glet and triplet spin configurations 0, [E IE] • State-of-the- 
art experimental techniques allow one to manipulate spin 
states of two electrons in a lateral, gate-defined DQD on 
time scales of hundreds of nanoseconds @, [1] , which has 
led to a renewed interest in factors that limit spin coher- 
ence. In theoretical investigations, the hyperfine interac- 



tion 



[E EE II > EH and phonon related phenomena me- 
diated by the spin-orbit coupling 
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are considered to be the major dephasing channels. This 
yields good agreement with experimental results for sin- 
gle electron spins, showing relaxation times of the order 



of milliseconds [13, 123, |21| ■ In a DQD system, the ex- 
perimental lower limit for coherence times lies in the mi- 
crosecond range This increased dephasing cannot be 
accounted for by the spin-related interactions mentioned 
above. 

In this Letter, we show that an efficient pure dephasing 
channel is always present in systems of two electron spins 
localized in coupled semiconductor QDs. This channel 
results solely from the charge-phonon interaction in the 
presence of inter-dot tunnel coupling, and is essentially 
due to the distinguishability of singlet and triplet states 
via Pauli-blocking of certain transitions in the triplet 
case. The key feature of this decoherence process is that 
it does not require any spin-environment interaction and 
relies only on the mechanisms (tunnel coupling and the 
Pauli principle) that are essential for the implementation 
of quantum gates. In particular, it appears also in mate- 
rials with negligible spin-orbit and hyperfine couplings. 

Qualitatively, in the lowest energy state of the two- 
electron system, each dot is occupied by a single elec- 



tron. The spin configuration of the system may then 
be either singlet or triplet. In the former case, the or- 
bital (spatial) wave function is symmetric and a transi- 
tion to a higher-energy doubly occupied state is possible. 
This is forbidden by Pauli exclusion in the triplet config- 
uration with an anti-symmetric orbital part. Although 
such transitions are completely inefficient at sub-Kelvin 
temperatures because of negligible occupation of the re- 
quired phonon states, a two-phonon process is still pos- 
sible, in which absorption of a p honon is followed by re- 
emission of another one 
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. In such a process, the 
doubly occupied state is involved only "virtually" and 
energy conservation requires only that the two phonons 
have the same energy (that is, the scattering is elastic) 
but this energy can be arbitrary. Therefore, even at low 
temperatures, phonons scatter on a DQD in the singlet 
state, while a DQD in the triplet state is transparent to 
phonons. This distinguishability leads to pure dephasing 
of any singlet-triplet superposition, in some sense anal- 
ogous to the "collisional decoherence" of the orbital de- 
grees of freedom [III, • Although the intensity of this 
process drops down at low temperatures because of de- 
creasing two-phonon spectral density at low frequencies, 
the temperature dependence is only polynomial (as op- 
posed to exponential suppression of real transitions). As 
we will show, at sub-Kelvin temperatures, at witch spin 
coherent control experiments on DQDs are performed, 
the two-phonon process can still lead to pure dephasing 
times as low as tens or hundreds of microseconds. 

In order to quantitatively estimate the effect of this 
dephasing we consider two electrons in laterally coupled 
quantum dots interacting with a phonon reservoir. The 
dots are considered identical and the model is restricted 
to the ground state in each dot. The Hamiltonian of the 
system is then 
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The first term describes the electrons and has the form 
Hbqd = ~h ^ ( a L a Rs + n - c - 



(1) 



ij.k.l 



where a,- s ,aj. are the electron annihilation and creation 

ta ' IS 

operators with i = L,R denoting the left and right dot, 
respectively, and s =Til labeling the spin orientation. 
The first term in Eq. JT]) accounts for single-particle inter- 
dot tunneling. The second term describes the Coulomb 
interaction, with V i]k i = Vjuk = Vidij = Vlkji (the 
wave functions may be chosen such that the matrix el- 
ements are real). For identical QDs the Coulomb ma- 
trix elements are also invariant under the interchange of 
the dots, L<->R. Among the Coulomb terms, Vlrrl = 
Vrllr = Ui and Vllll = Vrrrr = U 2 are the en- 
ergies of the singly- and doubly-charged configurations, 
Vlrlr = Vrlrl = E x are exchange energies, Vllrr = 
Vrrll = tc2 is the coupling between the doubly-charged 
configurations, while Vrlll and equivalent terms account 
for the coupling between the singly- and doubly-charged 
configurations and will be denoted by tc- 

The Hamiltonian of the phonon reservoir is given by 

Hph = J2k,\ ^fe,A 6 fe,A^fe,A) where bk,x,V kX are phonon 
annihilation and creation operators for a phonon from a 
branch A with a wave vector k and hojk t \ are the cor- 
responding energies. The electron-phonon interaction is 
described by 



fc.A 
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where F[f R (k) = F^(k) exp[±ik x D/2] are coupling 
constants and D is the inter-dot distance. We include the 
deformation potential and piezoelectric couplings. The 
coupling constants for the longitudinal and transverse 
acoustic phonon branches are, respectively 2(| 27|. 
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where A =l,ti ,ta refer to the longitudinal and two trans- 
verse acoustic phonon branches. Here e denotes the elec- 
tron charge, p c is the crystal density, v is the normaliza- 
tion volume for the phonon modes, d is the piezoelectric 
constant, £o is the vacuum permittivity, e s is the static 
relative dielectric constant and a is the deformation po- 
tential constant. The functions M\ depend on the ori- 

' E3- For the 



entation of the phonon wave vector [26 



zinc-blende structure they are given by 
M\(k) = 2 k x k y (e\,k)z + kyk z (e\ tk )x + k g k x (ex,h)y 



where k = k/k and e.\ & are unit polarization vectors. 
The form factors !F(k) depend on the wave function ge- 
ometry and are given by T(k) — J d 3 rip* (r)e lkr i(j(r), 
where ^>(r) is the envelope wave function of an electron 
centered at r = 0. 

We will use the basis composed of the three triplet 
states 
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and the three singlet states 
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The triplet states and the |( 1,1)5) singlet involve elec- 
trons occupying separate QDs and, therefore, have lower 
energies than the other two singlet states. 

The eigenstates of the Hamiltonian i?DQD are the three 
triplets with the energy U\ — Ex, the singlet |(— )5) with 
the energy Er_-\g = XJ<i — tci, and the two states 



\S+) = 
\S\ = 



1 
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[|(+)5)+e|(l,l)5)], 



](l,l)5)-£|(+)5)], 



where £ = 2 y/2t / (U + V U 2 + 8t 2 ) , with the cigcncncrgics 
E± = E±VU 2 + 8t 2 /2. RereE=(U2+t C 2 + Ui+Ex)/2, 
U = U 2 + t C 2 - Ut - E x , and t = V2(t c - h). In the 
weak tunneling regime, t <C U, one has 1 and |5+) f=a 
|(+)5), |5_) « |(1,1)5). The degenerate triplet states 
and the singlet state |5_) are the lowest energy states. 
In the following, phase decohercncc of a superposition 
of the |5_) singlet state and one of the triplet states is 
investigated. 

Since the electron-phonon interaction conserves spin, 
the singlet state |5_) is not coupled by phonon- assisted 
transitions to the triplet states. Calculation shows that 
|5_) is also decoupled from |5 + ), so the only nonzero 
off-diagonal matrix element of Hi nt involving 1 5_ ) is 



(5_|tf int |(-)5) = 
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For this coupling, the spectral density of the phonon 
reservoir (as defined, e.g., in 27J) takes the form 



the two relevant singlet states, 
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x [{n k ,\ + l)d(u - Wfc a) + n k 6(oj + 0Jk,\)] , 



where n k ,x is the Bosc distribution. 

The two electron system is descried by the reduced 
density matrix puqd = TrR p, where Ttr denotes the 
partial trace over the reservoir degrees of freedom and p is 
the density matrix of the complete system. Its evolution 
can be described using the time-convolutionless (TCL) 
projection operator method [28]. For factorized initial 
conditions (pure initial state), the TCL master equation 
takes the form d[Pp(t)]/dt = K.{t)Vp{t), where the pro- 
jection operator V is defined by Vp(t) = TrRp(i) ® pr, 
IC(t) is the TCL generator, and pr. is the density matrix 
of the phonon reservoir at the thermal equilibrium. Ex- 
panding the TCL generator up to the fourth order yields 
El £(*) = ZLiM*): with K x {t) = JC 3 (t) = and 



K 2 (t) = f dtiPC(t)£(Jti)P, 
Jo 
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P£(t)(£(ti)P£(t 2 )£(t 3 ) + C(t 2 )VC(h)C(t 3 ) 
£(t 3 )7>£(*i)£(* 2 ))7>] , (3b) 



where C(t) is the Liouville operator, C(t)p(t') = 

-(i/h)[H int (t), P (t')}. 

The coherence between the low-energy singlet and any 
of the triplet states is stored in the off-diagonal elements 
of the DQD density matrix, r, = (5_|/0dqd|(1, l)Ti) , 
i =T? 0. Upon explicitly evaluating the generators 
in Eqs. (f3"al b), the evolution equation for any of these 
singlet-triplet off-diagonal elements can be written in the 
form dri/dt — K(t)ri. Here, we are interested in the loss 
of coherence, that is, in the evolution of the modulus of 
This is given by 



d\rj_ 
dt 



= ReK{t)\n\. 



The function K(t) varies with time only for t < h/(k^T) 
(reservoir memory time), which is of the order of picosec- 
onds. Since the time scales relevant for the present dis- 
cussion are many orders of magnitude longer, we can 
safely use its long-time (Markovian) limit K = K(t = 
+oo). Then, K may be separated into two parts K — 
+ #(2), T he first term is 

Rc/i (1) = -tt(1 - e)R{-uj - AE), 

where hu>o = Er_\g — E- is the energy difference between 
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and the spectral density of the phonon reservoir R(uj) is 
given by Eq. This term describes the Fermi golden 
rule rate of real single-phonon transitions between the 
and |(— )S) states with fourth order corrections due 
to phonon-induced energy shifts and coupling renormal- 
izations. The second contribution is 



ReK^ 
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R(ujq — u>)R(—ojo +u>) + R(cvq + lo)R{—u)q — uj) 
~ 2 

R(ujq — U))R(— UJq + lo) — R(u>Q + Ll>)R(~Ll>o — uj) 
IjJoUJ 

and accounts for the two-phonon elastic scattering pro- 
cess. 

In the calculations, DQD geometry and material pa- 
rameters are taken which correspond to lateral, gate- 
defined QDs made in the two-dimensional electron gas 
(2DEG) of a doped GaAs/AlGaAs interface heterostruc- 
ture @, H3|. Two-dimensional Gaussian single electron 
wave functions are used with 170 nm full width at half 
maximum of the probability density. We set D = 200 nm, 
and U = 0.8 meV and use the material parameters 
q = 5100 m/s, c t = 2800 m/s, e s = 13.2, d = 0.16 C/m 2 , 
a = —0.8 eV, and p c — 5360 kg/m 3 . For simplicity, we 
use GaAs bulk phonon modes. 

The pure dephasing rates resulting from the two- 
phonon (scattering) process are shown in Fig. [TJa) as 
a function of temperature. Dephasing rates due to the 
single-phonon assisted transition for t = 0.3 meV are 
shown in the same figure for comparison. At low temper- 
atures, at which experiments on DQD ensembles are per- 
formed, the single-phonon transition is suppressed. At 
these temperatures, the dominating decoherence mecha- 
nism is the elastic scattering. This two-phonon process 
is much less influenced by decreasing the temperature 
since it involves only a virtual transition to a higher en- 
ergy (doubly charged) singlet state. The resulting pure 
dephasing rates at sub-Kelvin temperatures are relatively 
high, compared to experimentally achievable gate opera- 
tion times lZI . 
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FIG. 1: Two-phonon induced pure dephasing rates: (a) As 
a function of temperature for different tunneling parameters, 
(b) As a function of the tunneling parameter for different 
temperatures. Gray line in (a) shows the dephasing rates 
from single-phonon transitions for t = 0.3 meV. 



Fig. [T^b) shows the pure dephasing rates as a function 
of the coupling parameter t. This parameter, which is 
also crucial for unitary operations on the two-qubit sys- 
tem, affects the dephasing rate via the energy difference 
between the singlet states, Tiujq, and via the mixing pa- 
rameter £, which enters the spectral density in Eq. (J2]). 
Obviously, the dephasing vanishes for uncoupled dots. 
However, for non-zero coupling, the dephasing rate grows 
rapidly with t. 

In conclusion, we have shown that inter-dot tunneling 
and the Pauli principle, which are necessary for two-spin 
quantum gate operation, lead to singlet-triplet spin de- 
phasing in two-electron DQDs. This dephasing process 
does not depend on any interactions between spins and 
their environment and is therefore qualitatively different 
from the dephasing mechanisms discussed so far. Un- 
like single-phonon-assisted real transitions between sin- 
glet states, which are suppressed at low temperatures, the 
two-phonon elastic process remains non-negligible in the 
sub-Kelvin range relevant for spin coherent control exper- 
iments on DQDs. For a gate-defined GaAs/InGaAs DQD 
system, this process leads to dephasing rate of the order 
of 10 /is at T = 0.5 K. This is consistent with the rela- 
tively short dephasing times found in experiments 0, H[ • 
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